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We consider electron transport in a ballistic multi-mode 
channel structure in the presence of a transversely polarized 
far-infrared (FIR) field. The channel structure consists of 
a long resonance region connected to an adiabatic widening 
with a potential barrier at the end. At frequencies that match 
the mode energy separation in the resonance region we find 
distinct peaks in the photocurrent, caused by Rabi oscilla- 
tions in the mode population. For an experimental situation 
in which the width of the channel is tunable via gates, we 
propose a method for reconstructing the spectrum of propa- 
gating modes, without having to use a tunable FIR source. 
With this method the change in the spectrum as the gate 
voltage is varied can be monitored. 



I. INTRODUCTION 

Attempts have been made during the past few years 
to observe effects of the application of a far-infrared 
(FIR) field on electcoiis propagating coherently in narrow 
quantum channels au In focus of such investigations has 
been quantum point contacts (QPCs) defined by split- 
gate depletion of a two-dimensional electron gas in a 
GaAs/AlGaAs heterointerface.ErB So far the experiments 
employing both transverse and parallel polarization have 
been explained by heating and rectification effects, while 
clear evidence of FIR-ficld influence inside the channel is 
still missing. 

Considering the pronounced transverse quantization 
into waveguide modes it becomes particularly interesting 
to consider transverse polarization of the FIR field, since 
it allows strong mode coupling even for a homogeneous 
FIR field. For transverse polarization it is well known 
from theoretical considerations that transitions between 
modes take place mainly at points where the mode energy 
separation, E m — E n , equals the photon energy fiu nu At 
such points momentum conservation is possible. In the 
case of a QPC the mode energy separation changes and 
resonance conditions can be fulfilled only in a small re- 
gion, which reduces the sensitivity to the FIR field. This 
could explain why not even experiments using transverse 
polarization have succeeded in providing evidence of co- 
herent FIR-field influenced transport. 

The natural way to increase the sensitivity to a FIR 
field is to extend the region of resonance by considering 
a longer channel in which the width is constant. To con- 
sider long coherent channels has in the last few years hs=. 
come realistic due to progress in fabrication techniques. E£l 



Since FIR-ficld induced transitions between propagating 
modes alone will not change the transmission probability, 
it is necessary to incorporate some detection mechanism 
that discriminates between modes. This idea has been 
developed in the special case when only one mode enters 
the channel and a QPC is used for jdetection, which leads 
to a reflection of excited electrons OO 

In this work we consider the complementary case of 
a multi-mode channel. A new feature is that there is 
a possibility for cancellation if one electron climbs the 
mode spectrum while another one descends. Since there 
is no lower limit to the mode energy separation when 
the restriction of only one propagating mode is relaxed, 
a wider range of frequencies is of interest. For detection 
we use a combination of an adiabatic widening and a line- 
gate barrier, which leads to an enhanced transmission for 
excited electrons. We find a general expression for the 
scattering states, that includes all relevant mode mixing, 
by solving an eigenvalue problem. An explicit analytical 
expression is derived which is valid in the limit of weak 
FIR-fields, where resonances between pairs of transverse 
modes will dominate. We also propose an experimental 
method for finding the mode spectrum using only fixed 
frequency FIR-sources. 



II. MODEL SYSTEM 

We consider a quasi-onedimensional channel that 
smoothly connects two reservoirs. The channel consist 
of three regions, each with its own purpose (see Fig. [l]). 
First, there is a resonance region featuring strong size 
quantization, in which an external high frequency elec- 
tric field polarized across the channel can induce mode 
transitions. Second, there is an adiabatic widening in 
which there is a conversion of the energy stored in the 
transverse direction into kinetic energy for longitudinal 
motion. Finally, there is a barrier that blocks slow elec- 
trons. 

Experimentally such a system can be realized in differ- 
ent ways, but in order to see resonances fully developed 
it takes a system that preserves coherence all along the 
channel. We primarily have a split-gate channel in a 
GaAs/AlGaAs-interface in mind. The barrier can then 
be created by a line gate across the channel. Unlike a 
split-gate barrier it gives the same barrier for all modes. 
The height of this barrier is adjusted so that no transport 
takes place in the absence of the FIR-ficld. Then only 
those electrons that have been excited can pass. Since 
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there is mode pumping on one side only a photocurrent 
will be generated even in absence of a driving voltage. 
We are interested in this zero bias photocurrent. 



| E m (x) 




FIG. 1. Geometry of the channel and mode energies 
E m (x). Between x = and x = L there is a resonance region 
with an ^-independent mode energy spacing. At Xb there is a 
line-gate which gives rise to the barrier in the potentials. 



We assume that the width of the channel varies slowly 
on the scale of the electron wavelength which allows us. ta 
separate the transverse and the longitudinal motionJiMa 
Phase breaking processes and collective effects are ig- 
nored, and the current is calculated from single particle 
transmission probabilities in the spirit of the Landauer 
approach. E£l 



III. MULTI-MODE SCATTERING STATES 

We solve a single particle Schrodinger equation with 
the following time dependent Hamiltonian: 

H(t) = -^(p-^A(x,t)) 2 + U(x,y). (1) 

Here p is the momentum, m* is the effective mass of 
an electron in the two-dimensional electron gas and the 
potential U(x,y) confines the electron to the channel 
and reservoirs. The vector potential A(x, t) in Eq. (|l|) 
describes a homogeneous electromagnetic field of angu- 
lar frequency u> and amplitude E polarized in the in- 
direction: 



A(x, t) — — E cos(wi) e y . 



(2) 



We now make a separation Ansatz for the wave func- 
tion: 



(3) 



Here & m (x,y) are solutions to the transverse eigenvalue 
equation: 



h 2 d 2 

^7^2 $m ( a; ' V) + U ( x ' V)®rn(x, y) 
= E m {x)$ m (x,y) . 



(4) 



If the channel geometry changes slowly along the channel 
we can neglect ^-derivatives of <f> m (x, y) when insetting 
the Ansatz into the Schrodinger equation and find:E3 

Ti 2 d 2 d 

- ^r^-r*m(a;,*) + E m {x)m m (x,t) - m—^ m (x,t) 

2m ax A at 

= -2i cos(wi) ^2 M mm< (x)^m' t) , (5) 

m 

where the intermode transition elements M mm f(x) — 
-M m ' m (x) are given by: 



■M-mm' (#) 



TieE 
2m*u> 



f d 
J K l (x,y)^m'(x,y)dy . (6) 



We now look for solutions to Eq. (||) in the resonance 
region where U(x, y) is assumed to be x-independent and 
use the following Ansatz: 



i[Px-(E+lhu) t]/h 



(7) 



When this Ansatz is inserted into Eq. (||) we get an 
infinite set of coupled equations for the coefficients C m i : 

(P 2 /2m* -K ml {E))C ml (8) 

= — ^ iM mm r (C m /(J_1) + C, n '(/ + 1)) , 
n I 

where: 

K m i(E) = E + IKu - E m , (9) 
Eq. (|J) can be put in the form of an eigenvalue problem: 
Kmi,m'l' (E) C m , v = P 2 /2m*C ml , (10) 

where: 

K m i <m >i'{E) — K m i[E) 6 mm r6u> 

-iM mm ,(8 l{v+1) + <%/_!)) , (11) 

with 6 being the Kronecker delta. 

We will now use the fact that the resonance region 
has a finite length L, and that our solution inside the 
resonance region arise due to an incoming wave from the 
reservoir of the following form: 



* I n N ] f(x,y,t)=$ n (x,yy- 



J [J p n (x)dx-Etj /h 



VPn(x)/p n (0) 



(12) 



where p n (x) = yj2m*{E — E n (x)). This in combination 
with the restriction: \M\ <C Tiw, allows us to get a good 
approximation to the scattering state using a reduced 
number of states in the eigenvalue problem in Eq. ( |io| ) . 
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We pay attention to resonance effects only. For this 
reason we keep only one coefficient C m i for each mode, 
i.e. the one that gives the best kinetic energy match 
with the incoming state. Explicitly we define l mn as the 
integer that minimizes \E m — E n + l mn fiu>\, and keep only 
coefficients C m i mn . 

By this approach we miss the possibility of indirect 
resonant coupling via non-resonant states. On the other 
hand all ignored states necessarily give a kinetic energy 
mismatch of at least Tilu/2. This leads to a reduction of 
the effective coupling, via one or several such states, by 
at least a factor of the order \M\/%u). The characteristic 
length scale for population oscillation is at resonance: 
A_r = X-K/\M\, where Xis the de Broglie wavelength and 
K is the kinetic energy.EJ Obviously this length scale is a 
factor of %uj/\M\ larger for the indirect resonant coupling 
than for the direct one. Therefore there is a regime of 
weak FIR fields in which X 1 ^ 2> L, where it suffices to 
keep only the "most resonant" coefficients. 

We consider a situation in which the transport is adi- 
abatic outside the resonance region. This is realistic if 
the FIR field is sufficiently weak. Far from the resonance 
region, FIR field induced transitions are suppressed by 
a factor of the order \M\/hu>, because of a kinetic mis- 
match of the order Tllu. Near the resonance region the 
mode potentials are taken to be rapidly varying on the 
scale of Xr (but still slowly on the scale of A). In this 
case there is no room, for a resonance to develop a signif- 
icant change in the mode population. One can think of 
the FIR field as being suddenly switched on at x = and 
suddenly switched off at x = L. In particular we can as- 
sume that the incident mode, n, is still fully occupied at 
x = (where the resonance region begins) and that the 
mode population at x = L (where the resonance region 
ends) will remain for x > L. 

In the resonance region we get, using Eq. (||), the 
following expression for the scattering wave function: 



xe 



1 \ P n r ' E X-(E + hu>l nm )t \/h 



(13) 



where P 7 [% = yj2m*K% , and K% E and C^ m are 
the eigenvalues and the orthonormal eigenvectors of the 
following reduced eigenvalue equation: 



Y,[K m i mn {E)5n 



(14) 



(r) 

and the coefficients a n E are chosen such that only mode 
n is populated at x = 0: 



E n,E' 



(15) 



In the widening region [x > L) we get, using the as- 
sumption of adiabaticity, the following expression for the 
scattering wave function: 



m 

i\j Vn s E,rn{x)dx-{E+hojl mn )t\ /ft 

x$m(x,v)- / 7 w =7fr 

■\/Pn,E, m (.X)/p n>E>m {L) 



(16) 



where p n ,E, m (x) = ^2m*(E + hivl mn - E m (x)), and the 
transition amplitudes t mn are given by: 



iP {T> L/h 



(17) 



When we call t mn (E) transition amplitudes we implic- 
itly assume that scattering takes place between states 
with the same momentum. The coupling is strong only 
between states whose difference in kinetic energy AK 
fulfills AK < \M\. Thus for coupled states we have 
Ap/p w AK/2K < \M\/2K, and assuming |M| < K 
we can forget this difference in momentum. 

Our scattering solution \E' Tti £;(x, y, t) is valid for suffi- 
ciently weak FIR fields since it does not take reflections 
into account. In general there will be reflections both 
at the entrance and at the exit, but they will be sup- 
pressed if \M\ <§; K. In principle we could handle more 
than the weak FIR field limit by using an exact matching 
condition that takes into account also the derivative of 
the wave function, However, in the strong field limit one 
would also have to account for multiple reflection from 
the entrance, the exit and from the barrier. 

An electron near the band bottom will not fulfill the 
condition \M\ <C K. For such an electron we fail to make 
reliable predictions. On the other hand, since each energy 
interval is weighed equally in the Landauer formula we 
can ignore the contribution from slow electrons provided 
that the energy interval of integration is much larger than 
\M\. n 

According to the Landauer approachO we have, at zero 
temperature, the following expression for the current. 



2e 



T n {E) 



(18) 



For T n (E) which is the total transmission probability we 
have: 



T n (E) = W mn (E)\t mn (E)\ 2 e{E - E n ) 



(19) 



where W mn (E) is the barrier transmission probability for 
an electron in mode m at energy E + l mn huj. Assuming 
the barrier to be smooth we can ignore tunneling and use 
the following expression: 



W mn (E) =@[E + l nm htu - (E b + E m (x b ))} 



(20) 



Here E b is the potential energy due to the line-gate and 
E m {xb) is the residual effective mode potential induced 
by the transverse confinement in the widened region. We 
consider the case when E b is adjusted so that there is no 
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transmission in absence of the FIR field, and we ignore 
E m {xb). The part that is reflected from the barrier need 
not be accounted for since it will travel in the reverse 
direction until it reaches the reservoir. 



IV. PAIR-COUPLING APPROXIMATION 



We now restrict our attention to situations in which: 



\(E m , -E m )-(E m -E n )\ » |M| 



(21) 



this case, only pair coupling is important. If there are 
some modes: m', m and n that do not satisfy the condi- 
tion in Eq. (pl|), a pair-coupling approach is invalid for 
frequencies that couple these levels. 

Two different kinds of coupling can be distinguished 
in the pair-coupling approximation. There is an "ab- 
sorption coupling" in which an electron entering in mode 
n at an energy E is coupled strongly to a higher mode 
m at a higher energy E + huj. There is also an "emission 
coupling" in which the same electron is coupled strongly 
to a lower mode m at a lower energy E — Tilo. These cou- 
plings are not active simultaneously with our restriction 
in Eq. ( ^l|) on the spectrum. 

The "absorption coupling" and the "emission cou- 
pling" give the same eigenvalue equation by suitable sub- 
stitutions. For the "absorption coupling" we have l nn = 



K^ a (E) 



ABSi 
ran 

ABS 
ran 

ABS 



AK'„ 



(K ml 

K m l ~ 



-K nQ )/2 

K n0 



ABS 



K 



ABS 



(E), 



while for the "emission coupling" we have l r , 
Imn = — 1 and we introduce: 

K™{E) = (K m( _ x) +K n0 ) /2 



(22) 



and 



EM 



K™(E) . 



(23) 



Introducing M — M nm — —M rnn and dropping sub- 
scripts and superscripts for now, we get the following 
symmetric and _E-independent form of Eq. (^) : 



-AK/2 -%M \ f C„ 
iM AK/2 ) \ C„ 



K. 



(24) 



The eigenvalues are: = ± ^(AK/2) 2 + M 2 and the 
orthonormal eigenvectors are: 



(+) 



-ia 
b 



Cri 



(-) 



lb 

a 



(25) 



where: 



a = —= \J\ — sin a 
V2 

1 



b = —= v 1 + sin a 
V2 



sin a 



AK/2 



^{AK/2f + M 2 



(26) 



Using the eigenvectors in Eq. (|25|) we find the ex- 
pansion coefficients using the matching condition in Eq. 



,(+) 



ai ' = — ib 



(27) 



Then we have for the transition probability in Eq. (|l9|): 



\tmn{E)\ 



2 J 7mn Sin {(Jran (E)L) , K mn {E) > 

, K mn (E) < 

M 2 



(AK mn /2) z + Ml n 



Qmn(E) = 



P^e-pLe tt / (AK mn /2) 2 + Ml, 



A — Xmn(E) 



2h A V K*JE) 

h 



y/2m*K mn (E) 



(28) 



We have already assumed a large kinetic energy: 
K n o 3> |M|, in the incident mode. Then, if K < we 
know that the kinetic energy in the other mode must be 
negative. This implies a large kinetic energy mismatch: 
\AK\ ^> \M\, and thus corresponds to an off- resonant 
situation (7 <C 1). Thus we are allowed to ignore the 
coupling and put \t mn (E)\ 2 — when K < 0, instead of 
trying to find an exact expression. 

The approximation for q mn (E) in Eq. (|28| ) is good 
when (AK/2) 2 + M 2 < K 2 However, since K n0 > \M\ 
we find that the approximation for q m n(E) is poor only 
in off-resonant situations: (AK/2) 2 ^> M 2 . Therefore 
we can use this approximation generally. 

From Eq. ( |2"i| ) we see that the population oscillates 
between the incident mode n to the other mode m along 
the channel. This is similar to Rabi oscillation in time 
of a two-level system. The wavelength of oscillation is 
given by ir/q, and the resonance strength is given by 7. 
It is clear from the expression for 7 that the pair comes 
to resonance when: \AK\ < 2M . 

The following relations are easily verified: 



K 



ABS 



(E) 

^ABS 



AK Aaa = AK 

mn run 



K™(E + hu) 

EM 



which leads to: 

|trW = l«£ + Mi 2 



(29) 



(30) 



Equation ( |30| ) states that as the electron from \n, E) is 
pumped up there is another electron from |m, E + hio), 
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which is pumped down, completely canceling the net 
change in the mode population. This is provided that 
both the scattering states are occupied. The cancelation 
which is a manifestation of the orthogonality between 
scattering states, appears when all coupled modes are oc- 
cupied. For this reason only pumping between states that 
have different occupation factors will give a net oscilla- 
tion along the channel. At zero temperature this happens 
for electrons that are closer than fuj to the Fermi level 
(see Fig. |). 



Total energy 




Kinetic energy 



FIG. 2. In the limit of weak coupling and zero temper- 
ature, only pair transitions close to the Fermi level matters. 
In this case an analytical expression for the transition ampli- 
tudes can be found. 



We consider the case: T — 0, Eb = Ep and E m (xb) = 
0. In this case only electrons that have been excited into 
an empty state will be transmitted over the barrier. It 
is clear that only "absorption coupling" can lead to this 
and we find the following expression for the current: 



N M 



/ = EE J « 



(31) 



n— 1 rn > n 



with the partial currents given by: 



1 mn 



K* b (E)\*e(E-E n )dE, (32) 



Ep — huj 



In the sums in Eq. ( |3l| ) , N is the number of propagating 
modes and M can be chosen to fulfill: Em — En > hu>. 
For a symmetric confining potential there are no transi- 
tions for which m + n = even. The transition probabil- 
ities \t^ s {E)\ 2 are found from Eq. ( |28| ) where for the 
"absorption coupling" we have from Eqs. d2^) and (0) 



K mn {E)=E 



TlUJ - (E m + E n ) 



{Er, 



2 

En 



(33) 



From Eqs. (|3J), and (|33j) it is clear that the pair 
resonances give rise to peaks in the current for certain 
frequencies f mn — (E m — E n )/h. The width of these 



peaks is: A/„ 



2M mn /h. 



V. RESULTS 

We will here consider the case of a square well confin- 
ing potential with impenetrable walls at y — ±d/2. In 
this case we have from Eq. (||) the following transition 
elements: 



M mn = 



AfieE 



3mn 



Sujdm* 2(m 2 — n 2 ) 



Jm-\-n,odd 



First we will demonstrate the frequency and field 
strength dependence of the current found within the pair- 
coupling approximation (see Eq. (|3l|)). In Fig. (|) we 
plot the current for a channel for which L = 5 /im, N = 6 
and E F = 14meV. 

Photocurrent (uA) 

m 




Frequency (THz) 



FIG. 3. Plot of the current versus both v = u>/2ir and E 
for a 5 /im long channel with six propagating modes. Peaks 
(light areas) correspond to pair resonances. The right column 
gives the scale for the current in microamps. 



For weak FIR fields the pair-resonance peaks are well 
separated and the pair coupling approximation is good. 
However, in the upper-right part of the figure there is a 
slight overlap, and using the general scattering state in 
Eq. (|l3|) would lead to a more accurate result in this 
case. 

In the situation where the confining potential is tun- 
able by means of gates, it is interesting to consider the 
effect of varying the gate voltage, V g . The confining po- 
tential can change in different ways when V g is changed. 
As an example we will consider what happens if only the 
width d is changed, still using a square well potential. 
Then we have: 



E,, 



E f 



N 2 (V g ) ' 



(34) 



where N(V g ) still is the number of propagating modes. 
In order to make conncction-.to the famous split-gate ex- 
periment by van Wees et. aH, we use: N(V g ) = 20 + 91^. 

We can now observe the current peaks without sweep- 
ing the frequency. By instead sweeping the gate voltage 
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we change the mode energy separation E m — E n and thus 
the resonance frequencies. When E m — E n matches the 
applied frequency a pair will come to resonance. In Fig. |] 
we show the result of sweeping V g for a set of frequencies, 
when E = 5 V/cm. In this case one can verify that the 
condition on the spectrum in Eq. (21) is clearly fulfilled 
in the considered frequency range. 



AE (meV) 



2.0 



1.5 



1.0 



0.5 



hv=2.3 meV 



1 Li 



60 nA 



n=l n=2 n=3 n=4 



hv=1.6 meV 



hv=1.0 meV 



hv=0.6 meV 




-1. 



-1.6 -1.4 -1.2 



FIG. 4. By sweeping V g for a set of frequencies and in- 
terpolating between the positions of resonance peaks in the 
current one can find out how the mode energy separations 
AE n = E n+ i — E n change with V g . 



Each frequency gives its own set of peaks. It is clear 
that for a high frequency it takes a more negative V g 
to separate the mode energies sufficiently. Therefore the 
peaks are shifted to the left with increasing frequency in 
the diagrams. By drawing a line that interpolates this 
shift of the peaks we get a complete picture of how the 
mode energy separations AE n = E n+1 — E n changes with 
V g . By adding these separations we get the full mode 
spectrum relative to the lowest mode, E n (V g ) — Ei(V g ). 



VI. DISCUSSION 

In order to have a coherent influence on the propagat- 
ing electrons the coherence time in the far infrared source 
must exceed the passage time of electrons. A highly co- 
herent source is typically not widely tunable. In experi- 
ments ofr-the, considered kind, fixed frequency FIR lasers 
are usedJacm Therefore it is experimentally advantageous 
that we can reconstruct the mode spectrum using only 
fixed frequency sources, and also see how it changes with 
V g . By the interpolation method one can gain informa- 
tion about all occupied modes for a given V gi and not only 
about the one which is closest to pinch-off. Theoretical 



considerations suggest that there is a transition from a 
parabolic confinement towards a more square- well like, as 
the channel is made wider .tZI Such detailed information is 
not found by measuring only the conductance versus V g 
since, due to a variationJn the charge density, the spec- 
trum changes with V^Jj'tl For channels with a parabolic 
confining potential, V(y) = 77i*Wq2/ 2 /2, magnetotranport 



experiments can be used to determine ujqlj 

From a physical point of view our system is interest- 
ing in that it enables a directed acceleration of electrons 
caused by a high frequency field, via a two step pro- 
cess. First energy is absorbed in standing wave excita- 
tions across the channel. In this process a strong influ- 
ence is possible since both energy and momentum along 
the channel is conserved at resonance. Then the energy 
stored in the transverse direction is released in the widen- 
ing and adiabatically converted into kinetic energy along 
the channel. 

A nonzero temperature will destroy our predictions 
in two ways. First, the coherence length will be re- 
duced because of enhanced phonon scattering. In exper- 
iments clear quantized conductance in a 5 /im channel 
has been demonstrated at 1.3 K.t2l Since such temper- 
atures are used also in FIR-transport experiments this 
seems within reach.B The second degradation comes from 
thermal smearing. As a criterion we can take that the 
temperature must be much smaller then the inter mode 
separation, which is typically of the order 20 K. 

We have ignored collective effects. In absorption expex, 
iments on wire arrays one finds a depolarization effect S3 
As a result the mode spectrum found from absorption 
spectroscopy can differ from that found from magnetore- 
sistance oscillations Ej This is not immediately general- 
ized to a non-homogeneous system. A difference is that 
we consider coherent transport and that the inhomogene- 
ity acts as a boundary condition for propagating elec- 
trons. We can not rule out a significant influence due to 
depolarization, but leave it for future investigations. 

Another collective effect is the static spatial modula- 
tion of the potential along the channel because of changes 
in the mode population. However such an influence is 
suppressed in a multi mode channel since only a frac- 
tion of the electrons are resonantly coupled and elec- 
trons moving with different velocities give depopulation 
at different locations along the channel. In addition most 
charge comes from slow electrons, and for these strong 
cancellation can be expected because of a strong energy 
dependence. 

In conclusion we consider the influence of a FIR field on 
coherent transport in a multi mode channel. We focus on 
a long channel in which resonance conditions prevail over 
a considerable distance in order to achieve a high sensi- 
tivity to an applied FIR field. We derive both a general 
scattering state and an approximation valid in the limit 
of weak fields, where pair resonance will dominate. In 
this limit the current will peak when Tilo = E m — E n . We 
also propose a way to analyze the spectrum, using fixed 
frequencies only. 
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